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We derive an energy-dependent decay-time distribution function from the multi-exponential decay of the ensemble 
photoluminescence (PL) of InxGai.xN/GaN quantum dots (QDs), which agrees well with recently published single- 
QD time-resolved PL measurements. Using eight-band kp modelling, we show that the built-in piezo- and pyroelectric 
fields within the QDs cause a sensitive dependence of the radiative lifetimes on the exact QD geometry and composi- 
tion. Moreover, the radiative lifetimes also depend heavily on the composition of the direct surrounding of the QDs. 
A broad lifetime distribution occurs even for moderate variations of the QD structure. Thus, for unscreened fields a 
multi-exponential decay of the ensemble PL is generally expected in this material system. 
(accepted at Physica Status Solidi (b). ©2008 WILEY- VCH) 



1 Introduction The InxGai.xN/GaN system has evolved 
into one of the most important material systems for solid- 
state light emitters. Applications include green and blue 
light emitting diodes (LEDs) |[l]|2l, laser diodes (LDs) Gl, 
and white light emitters ||4]|5|. Already at early stage of 
research the enormous efficiency of InxGa^xN optoelec- 
tronic devices has been attributed to zero-dimensional lo- 
calization centers [quantum dots (QDs)] within the In^GaixN 
layer(s) |l6]|7l. Despite tremendous advances in this research 
field, the understanding of essential properties of this ma- 
terial system still needs fundamental improvement in order 
to further advance InxGai_xN/GaN-based optoelectronic de- 
vices. In particular, the influence of the large built-in piezo- 
and pyroelectric fields in InxGai^N-based heterostructures 
on the emission processes is still not understood in de- 
tail. Interestingly, in InxGai_xN QDs the quantum-confined 
Stark effect (QCSE) and the corresponding blue-shift at 
high charge-carrier concentrations has been reported to be 
efficiently suppressed JS]. Still, the origin of the InxGai.xN- 
typical multi-exponential photoluminescence (PL) decay is 
a matter of active debate ll9l [T0lfTnfT2l . It has been pointed 
out by Bartel et al. lfT2l that the multi-exponential decay 
can be consistently explained with broad distribution of 
single-QD decay times. In this paper we show that such 
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a broad lifetime distribution — and, consequently, a multi- 
exponential decay of the ensemble PL — is the effect of the 
built-in piezo- and pyroelectric fields within the QDs and 
can be generally expected in this material system. 

2 Sample Preparation and Experimental Results 

The samples investigated in this work have been grown on 
Si(lll) substrate by low-pressure metal-organic chemical 
vapor deposition (MOC VD) using a horizontal AIX200 RF 
reactor. An AlAs layer was grown and subsequently con- 
verted to AIN as a nucleation surface |13||. In the follow- 
ing step an Alo.osGao.gsN/GaN buffer layer was grown at 
T=1150 °C up to a total thickness of 1 /im. The InxGai^N 
layer was grown at 800 °C with a nominal thickness of 
2 nm using trimethylgallium, trimethylindium, trimethyla- 
luminun, and ammonia as precursors. The QDs are formed 
by alloy fluctuations within the InxGai-xN-layer 1141 . The 
growth was finished with a 20 nm GaN cap layer grown 
during a heat-up phase to 1 100 °C. The PL was excited at 
353 nm by the second harmonic of a mode-locked Ti:sapphire 
laser. The temporal width of the laser pulses was 2 ps at 
a repetition rate of 80 MHz. The measurements were per- 
formed in a helium-flow microscope cryostat at a tempera- 
ture of 4 K. The luminescence was collected through a mi- 
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Figure 1 /^-PL spectra of a single InxGai^xN/GaN QD. 
(a) Time-integrated spectrum, (b) Time-resolved PL of the 
emission line indicated by the black arrow in (a). The 
decay is well described mono-exponentially (solid black 
line). Ten different single-QD transients have been inves- 
tigated by Bartel et. aJ lfT2l . all of which showed a mono- 
exponential decay. 
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Figure 2 Excitation-density dependence of the TRPL sig- 
nal at fixed detection energy (3.05 eV; black arrow in the 
inset). A stretched-exponential function with identical pa- 
rameters describes all transients (white lines). The inset 
shows the time-integrated PL of the sample. 



croscope objective. The detection system consisted of two 
0.35 m McPherson monochromators in subtractive mode 
and an ultrafast photo detector (multi-channel plate) pro- 
viding a spectral resolution of about 0.6 me V and a tem- 
poral resolution better than 30 ps. The single-QD measure- 
ments are described in detail elsewhere lfT2l . The PL spec- 
trum of the entire QD ensemble has its maximum at 3.05 eV 
and a full width at half maximum of « 75 meV (inset in 
Fig. |2]i. Sharp emission lines of discrete QD states were 
observed from 2.8 to 3.2 eV proving the QD origin of the 
entire emission 01211141 . 

Photoluminescence transients measure on emission lines 
from single QDs mono-exponential (Fig.[T]i ||T2l. The ensemble 
PL decay, in contrast, is multi-exponential for all detection 
energies. A similar behaviour — multi-exponential decay of 
the macro PL and mono-exponential decay in single-QD 
measurements — has also been reported by Robinson et al. 
ifTsl . Rice et al. {TSj have even convincingly demonstrated 
that they observe the decay of single excitons and biexci- 
tons in their experiments on single InxGa^xN QDs. More- 
over, the shapes of the ensemble-PL transients do not de- 
pend on the excitation density (Fig. |2|i, indicating that the 
multi-exponential decay can be attributed to a broad distri- 
bution of decay times lil2J rather than dynamical screening 
effects Ellis]. 

3 Distribution of Excitonic Lifetimes All QDs with 
the same excitonic transition energy E can be considered 
a subensemble within the entire QD ensemble. The distri- 
bution of decay times r within each subensemble shall be 
described by the distribution function /^(t). This distri- 



bution function can be approximated from TRPL measure- 
ments of the subensembles: To eliminate the noise of the 
experimental data, the multi-exponential subensemble PL 
decays are approximated by the Kohlrausch (or stretched- 
exponential) function |[T9l (solid white lines in Fig.|2]i 



lE{t)=lE,0CXp[-{t/T*EY 



(1) 



Here, t|; and (3e are the time and stretching parameters for 
the given energy E. t'J^ is identical to the decay time of the 
system te only if Pe = 1. It has no direct physical mean- 
ing for (3e j^ I. The decay of the PL signal at eleven dif- 
ferent equidistant detection energies between 2.82 eV and 
'3.18eV has been analyzed. The shapes of the transients 
depend on the detection energy, resulting in varying fit pa- 
rameters for different QD subensembles (0.35 ns< t^ < 
0.85 ns; 0.55 < /3e < 0.84). 

The PL decay of each QD subensemble can also be 
expressed as an integral over the (exponential) single-QD 
PL decays of all QD which form the subensemble. Thus, 
using the Kohlrausch function to describe the decay, we 
can expand Eq.[T]to 



lE{t) oc exp[-(i/T^)''^] oc / /b(t) exp{-t/T)dT 



(2) 
Mathematically, the subensemble decay Ie (t) is, hence, 
the Laplace transform of the lifetime distribution function 
fE{T) within the subensemble, which can, in turn, be ob- 
tained by an inverse Laplace transformation of Ie (t) 

m. 
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Figure 3 (a) Lifetime distribution functions /^(t), ob- 
tained by inverse Laplace transformation of the multi- 
exponential PL decay of different InGaN-QD subensem- 
bles (different detection energies E). The integral of each 
/£;(t) is normalized to I. (b) f{E, r) as gray-scale plot to- 
gether with the decay times of single InGaN QDs reported 
in Ref. lfT2l . The single-dot lifetimes agree well with the 
obtained distribution function. 



The results of the inverse Laplace transformation are 
shown in Fig. [3] A broad distribution of lifetimes for all 
transition energies is found, in good agreement with the 
single-QD TRPL results from Ref. |12] [black symbols in 
Fig. [2b)], which also show an appreciable scatter for iden- 
tical transition energies. All single-QD time constants fall 
in the range covered by the lifetime distribution function. 

For GaN/AlN QDs a drastic increase (several orders of 
magnitude) of the excitonic lifetimes for lower transition 
energies has been reported il22l . which is caused by the in- 
creasing QCSE for increasing QD height. Such an effect 
can not be observed in the lifetime distribution function in 
Fig. [3] The maximum of the distribution function is con- 
stant at « 1.2ns up to transition energies of w 3.05eV;for 
higher transition energies it shifts slightly to shorter life- 
times, accompanied by a broadening of the distribution. As 



will be alluded to in Sec. |4] the broadening of the distri- 
bution at higher transition energies can be explained with 
the delocalization of charge carriers that are only weakly 
bound to shallower localization centers. 

4 Quantum Dot Structure and Radiative Lifetimes 

The broad lifetime distribution within the InGaN-QD en- 
semble can be understood in terms of varying electron-hole 
wave-function overlap in different QDs. We will show here 
that the variation of the overlap is caused by differences in 
the built-in piezo- and pyroelectric fields within each local- 
ization center 

We have calculated the radiative lifetimes of different 
InGaN QDs using a three-dimensional eight-band k • p model. 
The model accounts for strain effects, piezoelectric and 
pyroelectric polarization, spin-orbit and crystal-field split- 
ting, and coupling between the valence bands (VBs) and 
the conduction band (CB). Excitonic corrections have been 
included using a self-consistent Hartree (mean field) scheme. 
A detailed description of the method can be found else- 
where |f23l. The radiative lifetimes Trad of the confined ex- 
citons have been calculated by ll24l 



3t;2 



Trad 



27reor7T.oCg?i 



(3) 



Here, Eq is the permittivity of free space, mo the free- 
electron mass, Co the vacuum speed of light, h the reduced 
Planck constant, and E^^ the transition energy of the exci- 
ton. n is the refractive index of the matrix material (GaN), 
which can be described as a function of the emission wave- 
length by a Sellmeier-type law 1251 . The effective oscil- 
lator strength /eft has been calculated by integrating the 
(anisotropic) oscillator strength /e over the unit sphere: 



^^«=4^ 



/ede 



/e 



2r 



dO 



mo Eg 



m\e-p\m' 



(4) 

Here, |ifj) are the kp electron and hole wave functions, 
which consists of a sum over all k-p basis states \i) (Bloch 
functions) multiplied with the respective envelope func- 
tions ipi: |ifj) = J2i=\ "^iN)- When calculating the matrix 
elements contribution arising from the Bloch parts and the 
envelope-function parts of the wave functions have been 
accounted for e is a unit vector indicating the polarization 
of the light. 

In agreement with our recent publications II26II27I . the 
QDs have been modeled as ellipsoids with height h and 
lateral diameter d. They are embedded in a 2 nm thick In- 
GaN quantum well (QW) with indium concentration a;„. 
The indium concentration within the QDs increases lin- 
early from the indium fraction of the surrounding QW (x„) 
to the maximum indium concentration x^ at the QD cen- 
ter Starting with a QD with a height of ft, = 2 nm, a lat- 
eral diameter of d = 5.2 nm, and indium concentration of 
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Xc = 0.5, the influences of three different structural param- 
eters on the radiative excitonic Hfetimes have been inves- 
tigated: The QD height (with h varying between 1.2 and 
2.8 nm), the lateral diameter (with d varying between 2.8 
and 7.6 nm), and the indium concentration within the QD 
[Xc between 0.3 and 0.6). All three series have been calcu- 
lated with the QDs embedded in an Ing iGao.gN QW, in an 
Ino.05Gao.95N QW, and directly in the GaN matrix, with- 
out a QW. In the latter case QDs with slightly higher In 
concentration (x^ ~ 0.7) or slightly larger diameter {d ~ 
8.8 nm), respectively, have been included in order to cover 
the experimentally observed energy range. Note that, due 
to the concentration gradient inside the QDs, the average 
indium concentration within the QDs .Xavg is much lower 
than the maximum indium concentration of Xc ~ 0.3-0.6 
(0.7) at the QD center A maximum indium concentration 
of Xc = 0.3 corresponds to an average concentration of 
Xavg ~ 0.08-0.15 depending on the indium concentration 
in the QW (x„ = 0.0-0.1). An x^ of 0.6 (0.7) corresponds 
to lavg « 0.15-0.23 (0.18-0.25). The magnitude of the de- 
viation Ax = Xavg — Xw ~ 0.08-0.15 is in good agreement 
with typical values for alloy fluctuations in InGaN QWs 



Neglecting the built-in piezo- and pyroelectric fields 
[Fig. 3(a)] all QDs show similar excitonic lifetimes around 
0.8-1. ns. The radiative lifetimes are insensitive to the ex- 
act QD shape and composition. Such a lifetime distribution 
would lead to mono-exponential PL decay for all transition 
energies. Indeed, such a behaviour has been observed by 
Senes et al. 1291 for In^Ga^xN-QDs grown by molecular 
beam epitaxy (MBE), were field effects are reported to be 
neghgible |[29| . 

For the MOCVD-grown QDs investigated in this work 
a strong deviation from the mono-exponential decay has 
been observed for all transition energies. This deviation 
can be explained if the built-in electrostatic fields are in- 
cluded in the calculations [Fig. 3(b-d)]. The radiative ex- 
citonic lifetimes scatter across a wide range between 1.0 
and 5.5 ns. Pronounced dependencies on all structural pa- 
rameters can be observed. Interestingly, the radiative life- 
times significantly decrease, when the indium concentra- 
tion of the surrounding QW is decreased from 10% to 
5 %, or even to %. Then, the matrix exerts a increasing 
compressive strain to the z-axis ([0001] -axis) of the QD 
which reduces the built-in field from « 2.5 — 4.3MV/cm 
if embedded in a 10 % QW down to w 1.7 - 3.1 MV/cm 
(« 1.0 - 2.6 MV/cm) in a 5 % (0 %) QW. The separation 
of the centers of mass of the electron and hole wave func- 
tions decreases from « 0.9 — 1.6 nm to ?a 0.7 — 1.2 nm 
(w 0.4 — 0.8 nm). Thus, not only the exact geometry and 
chemical composition of the QDs themselves, but also the 
properties of the direct surrounding area of the QDs are de- 
cisive for the radiative lifetimes. A large variety of different 
time constants can be found even for identical transition 
energies. For instance, among all considered model QDs, 
six different QDs emit at « 3.1 cV (Fig. 3). Although all 
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Figure 4 Radiative excitonic lifetimes vs. transition ener- 
gies calculated with eight-band k • p theory: (a) All model 
QDs considered in this paper (see text), neglecting the 
built-in piezo- and pyroelectric fields, (b-d) As a function 
of the QD height, diameter, and composition, including 
piezo- and pyroelectric effects. 



six QDs have nearly the same transition energy, their radia- 
tive lifetimes scatter appreciably between 1.3 ns and 4.2 ns. 
Thus, luminescence at a certain detection energy originates 
from a subensemble of QDs, all of which have the same 
transition energy, but significantly different excitonic life- 
times. 

The calculated lifetimes describes the experimental scat- 
ter well qualitatively, but quantitatively they are generally 
larger than the measured ones. On the one hand the un- 
certainties of the material parameters of InxGai xN pro- 
vide a reasonable explanation for this systematic devia- 
tion. On the other hand, Narvaez et al. ||30l have shown 
for InGaAs/GaAs QDs that the lifetimes of charged ex- 
citons (positive or negative trion) are shorter than that of 
the exciton by a factor of « 0.5. Therefore, another possi- 
ble explanation is that the PL originates from the decay of 
charged excitons rather than neutral ones. 

5 Conclusion We have extracted the photon-energy- 
dependent decay-time distribution function from the multi- 
exponential decay of the ensemble PL of In^Gai.xN/GaN 
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QDs, which agrees well with recently publish single-QD 
decay times. We have calculated the radiative lifetimes of 
localized excitons in InxGai^N/GaN QDs and shown that 
the built-in piezo- and pyroelectric fields within the QDs 
are the origin of the broad lifetime distribution: They cause 
a sensitive dependence of the radiative lifetimes on the QD 
geometry and composition. The lifetimes are also very sen- 
sitive to the chemical composition in the direct vicinity of 
the dots. Therefore, a broad distribution of excitonic life- 
times and, consequently, a multi-exponential decay of the 
ensemble PL is generally expected in this material system 
in the case of unscreened fields. A mono-exponential decay 
of the ensemble luminescence, on the other hand, indicates 
vanishing field effects. 
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